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Abstract. 

Context. Gas ionized by starburst activity radiates radio recombination lines (RRLs), from which one can derive 
the plasma conditions and the number of massive stars formed in the burst, free of dust obscuration effects. 
Aims. We aimed to find detectable RRL emission from additional extragalactic starburst systems and to use the 
line properties to estimate the properties of the ionized gas. 

Methods. We conducted a search for RRLs in the nearby extragalactic starburst or Seyfert galaxies NGC 1808, the 
Circinus galaxy, NGC 4038/9, IIZw40, NGC 6221, NGC 7552, IRAS 18325-5926, IC5063, and VV114. We used 
the Very Large Array with resolution of 3" to 32" and the Australia Telescope Compact Array with resolution 
of 10" to search for the RRLs H91q and H92a with rest frequencies of 8.6 GHz and 8.3 GHz. From the new 
detections we derive conditions in the starburst regions. 

Results. We detected for the first time RRLs from the starburst nuclei in the Circinus galaxy and NGC 1808. 
The Circinus galaxy was detected in RRL emission with a line strength integrated over the source of 3.2 mjy, 
making it the fourth-strongest extragalactic RRL emitter known at this frequency (after NGC 4945, M82, and 
NGC 253) and so is suitable for detailed study. The line and continuum emission from the Circinus galaxy can be 
matched by a model consisting of a collection of 50 to 10 000 H II regions with temperatures of 5000 K, densities 
of 500 cm~^ to 50 000 cm"'', and a total effective diameter of 3pc to 50 pc. The Lyman continuum production 
rate required to maintain the ionization is 1 x 10^^ s"^ to 3 X lO''^ s-\ which requires 300 to 9000 05 stars to be 
produced in the starburst, inferring a star formation rate of 0.2 yr~' to 6 Mq yr~'. NGC 1808 was detected 
in RRL emission at 3.9o" with a line strength of 0.47 mJy. No radio recombination lines were detected from the 
other galaxies surveyed to a 3cr limit of 0.3 mJy to 1.4 mJy. 

Conclusion. We have detected RRLs from two galaxies, adding to the small but growing number of known extra- 
galactic RRL emitters. The Circinus galaxy is strong and especially suited to high-quality follow-up spectroscopic 
study. We derived conditions and star formation rates in the starbursting regions. Uncertainties can be reduced 
by future multi-transition studies. 

Key words, galaxies: individual: NGC 1808 - galaxies; individual: Circinus Galaxy - galaxies: nuclei - radio lines: 
galaxies 



1. Introduction 

Starburst activity has a profound impact on galaxies, 
through prodigious formation of stars and clusters, ioniza- 
tion of gas and outflows driven by overpressure in the star- 
bursting region that redistribute the interstellar medium 
(ISM). Studies of the gas ionized in the burst are interest- 
ing because they provide information on the ionizing pho- 
ton production rate and hence the number of high-mass 



stars formed and on the conditions in the surrounding 
ISM. 

However, optical and infrared diagnostics of the ISM 
are hampered by dust obscuration associated with the 
molecular clouds from which the stars form. Radio stud- 
ies are free of dust obscuration effects and the detection 
of radio recombination lines is particularly useful as their 
strengths provide diagnostics of the plasma conditions and 
give dynamical information with arcsec resolution. 



* deceased 



RRLs were first detected from the starbursts in M82 
and NGC 253 by Shaver et al. (1977) and Seaquist & Bell 
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(1977) soon after the potential to detect and interpret 
extragalactic RRLs was shown by Shaver (1978). Those 
galaxies have since been stiidicd over a wide range of 
frequencies, giving constraints on the physical state and 
kinematics in the nuclear regions (e.g. Anantharamaiah & 
Goss 1997; Rodn'gucz-Rico ct al. 2006). 

However, surveys to find RRLs in other galaxies pro- 
duced no further detections for a period (Churchwell & 
Shaver 1979; Bell & Seaquist 1978; Bell ct al. 1984). When 
survey sensitivities were improved by an order of magni- 
tude using the Very Large Array (VLA) during a renewed 
effort in the early 1990s, RRLs were detected near 8.6 GHz 
from several bright starburst galaxies. Those new de- 
tections are NGC 660, NGC1365, NGC 2146, NGC 3628, 
NGC 3690, NGC 5253, M83, IC694, Arp220, Hcnize2-10 
(Anantharamaiah et al. 1993; Zhao et al. 1996; Phookun 
et al. 1998; Mohan et al. 2001), and NGC 4945 at mm 
wavelengths (Viallefond, private communication). 

During an extension of the RRL surveys to southerly 
starburst and Seyfert galaxies using the Australia 
Telescope Compact Array (ATCA), wc made three new 
detections; NGC 3256, the Circinus galaxy, and NGC 4945. 
The detection of NGC 3256 was reported in a previous 
paper (Roy et al. 2005a). Here, we report the detection 
of the Circinus galaxy and four upper limits (NGC 6221, 
NGC 7552, IC5063, and IRAS 18325-5926). The detection 
of NGC 4945 will be reported in a later paper. A short 
report of all three detections was published by Roy et al. 
(2005b). 

During an extension of the RRL surveys using the 
VLA, we made a new detection of NGC 1808 (Mohan 
2002) and established sensitive upper limits on four other 
galaxies (NGC 4038/9, NGC 7552, nZw40, and VV114), 
all of which are reported here. 

We adopt = 75 km s~^Mpc~^, qo = 0.5, and A = 
0, and give velocities in the heliocentric frame using the 
optical velocity definition throughout this paper. 



Haddad 1992) in which the continuum is determined for 
each baseline by a linear fit to the visibility spectrum. The 
final continuum and line images were made using natural 
weighting of the {u, v) data to achieve maximum possi- 
ble signal-to-noise ratio and we averaged together the two 
transitions and two polarizations. 

Uncertainties on the absolute flux densities have 
an 11% rms random multiplicative component due to 
flux-density bootstrapping and atmospheric opacity, a 
0.12 mjy rms random additive component due to thermal 
noise in a 1 MHz channel and a systematic multiplicative 
component of 11 % rms due mainly to the uncertainty in 
the Baars et al. flux-density scale. 

2.2. VLA Observations 

NGC 1808, NGC 4038/9, nZw40, NGC 7552, and W 114 

were observed with the Very Large Array (VLA) . The ob- 
serving parameters and results are summarized in Table 
2. Data reduction was similar to the ATCA reduction. 

For NGC 1808, the CnB data suffered from a high 
residual bandpass calibration error causing low SNR spec- 
tra. Thus, we based the analysis on the DnC spectral data, 
but show the continuum image from the higher resolution 
CnB data. 



2. Observations 

2.1. ATCA Observations 

We observed the Circinus galaxy, IC 5063, IRAS 18325- 
5926, NGC 7552, and NGC 6221 with the Australia 
Telescope Compact Array (ATCA). The observing param- 
eters and results are summarized in Table 1. 

Calibration and imaging were done using the AIPS 
software using standard methods. The flux-density scale 
assumed that PKS B1934-638 had a flux density of 2.99 Jy 
at 8295 MHz and 2.87 Jy at 8570 MHz, based on the Baars 
et al. (1977) flux-density scale. A phase calibrator was ob- 
served every half hour to correct the instrumental phase 
response. A bandpass calibrator was observed every few 
hours for correcting the instrumental frequency response 
(bandpass). Phase corrections obtained from self calibra- 
tion of the continuum source were applied to the spectral 
line data. Continuum emission was subtracted from the 
line data using the method UVLSF (Cornwell, Uson & 



Table 1. ATCA Observational Parameters and Results. 





Circinus galaxy 


NGC 6221 


NGC 7552 


IC 5063 


IRAS 18325-5926 


Array configuration 


750D (1993jul),750C (1994oct) 


750C 


750C 


750A 


750A 


iJats oi observation 


Ayyojui, lyy^tocr^u o6 nov^o 


lyy^octzz, 


lyy^octzz, zo 


lyy^mar.}, 4, o, o, * 


lyy^mard, 4, o, o, i 


Integration t iinc { ii } 


37 


8 5 


6 2 


9.7 


10 6 


Tiransitions 


H92q & H91q 


H92q & H91q 


H92q & H91q 


H92q &: H91q 


H92a & H91a 


i^rcst of RRL (MHz) 


8309.38 & 8584.82 


8309.38 & 8584.82 


8309.38 & 8584.82 


8309.38 & 8584.82 


8309.38 & 8584.82 


l^band-ccntrc (MHz) 


8295.00 & 8570.00 


8265.00 & 8540.00 


8265.00 & 8540.00 


8225.00 & 8487.00 


8141.00 & 8411.00 


Kystemic,helio, optical (kmS^"^) 

Distance (Mpc) 


480 


1477 


1590 


3402 


6065 


4.2 


18 


21 


45 


81 


Bandwidth (MHz), channels, IFs 


64, 64, 4 


64, 64, 4 


64, 64, 4 


64, 64, 4 


64, 64, 4 


Spectral resolution (kms~^) 


36 


36 


36 


36 


37 


Polarization 


dual linear 


dual linear 


dual linear 


dual linear 


dual linear 


Minimum baseline (kA) 


0.86 


1.27 


1.27 


2.13 


2.13 


Beam size (natural weight) ("x") 


11. 0" X 9.4" 


11.8" X 7.9" 


14.0" X 8.3" 


6.4" X 5.8" 


6.3" X 5.6" 


Beam position angle 


-76° 


32° 


16° 


24° 


.50° 


PViatjia (-aliVirafnr f'Rl Qi^n^ 




1740-.517 


2313-438 






Bandpass calibrator (B1950) 


0537-441, 1921-293 


1921-293 


0316-1-413, 0537-441, 
1921-293 


0537-441, 2251 + 158 


2251 + 158 


Line and Continuum Properties 












Peak continuum surface brightness (mjybeam"^) 


98.3 ± 15 


20.2 ± 2.3 


44.7 ±5.0 


180 ± 20 


30.9 ±3.5 


Total continuum flux density (mJy) 


234 ± 37 


20 ±3 


57 ±6 


177 ± 20 


30 ±4 


Noise per image channel (mjybeam"^) 


0.12 


0.35 


0.45 


0.33 


0.27 


Integrated line flux density (mJy) 


3.2 ± 0.6 


< 1.0 


< 1.4 


< 1.0 


< 0.8 


Integrated line flux or 3(7 upper limit (Wm~^) 


(20.1 ± 3.1) X 10"^^ 


< 1.0 X 10-^^ 


< 1.4 X IQ-^'* 


< 1.0 X 10-23 


< 0.8 X 10-23 


Measured line width (FWHM; kms""'") 


278 ± 10 










Deconvolved line width (FWHM; kms~^) 


260 ± 10 










Centroid helio. optical vel. (kms~^) 


464 ± 10 










Cont. flux density over line region (mJy) 


225 ± 35 










No. of beam areas where line is observed 


1.9 










Total 2.7 GHz continuum (mJy) 




220'' 


220" 


800° 


115° (at 2.4 GHz) 


Total 5.0 GHz continuum (mJy) 




110"* 


88^ (at 4.8 GHz) 


420" 





" Wright & Otrupcek (1990) 

Forbes et al. (1994b) 
° Roy (1995) 

Wiiitooak (1970) 
' Wrigiit (1974) 

^ Upper limits are three times the rms noise in a single channel. 

^ Upper limits assume a FWHM equal to the velocity resolution of the data. 



Table 2. VLA Observational Parameters and Results 



NGC 1808 NGC 4038/9 II Zw 40 NGC 7552 VV 114 



Arr&y configur&tion 


CnB 


DnC 


C 


DnC 




CnB 


D 


j5 


D&tc of obscrv&tioH 


1998novl0 


1999febl8 


1998dec27 


1999feb20 


2000augl4, sep29 


1998novl6 


1999apr23 


1999apr03 


Integration time (h) 


6.3 


6.3 














Transitions 


H92a 


H92a 


H92a 


H92a 


H92a & H93a 


H92a 


H91a 


H92a 


i^rest of RRL (MHz) 


8309.38 


8309.38 


8309.38 


8309.38 


8309.38, 8045.60 


8309.38 


8584.82 


8309.38 


I'band-centre (MHz) 


8285.100 


8285.100 


8264.725 


8264.725 


8288.185, 8025.080 


8264.900 


8535.100 


8288.185 


Kiyatomicholio, optical (klHS^^) 




1014 




1663 


789 


1653 


6016 


Distance (Mpc) 




11 




19.2 


10.5 




21 


77 


Bandwidth (MHz), diannels, IF 


46.9, 15, 1 


46.9, 15, 1 


23.44, 15, 2 


23.44, 15, 2 


23.4, 15, 2 


46.9, 15, 1 


46.9, 15, 1 


23.4, 15, 2 


Spectral resolution'* (kins~^) 


226 


226 


113.4 


113.4 


113 


226 


226 


115 


Polarization 


dual circular 


dual circular 


dual circular dual circular 


dual circular 


dual circular 


dual circular 


dual circular 


Minimum baseline (kA) 


3.8 


1.2 


2.0 


1.3 


0.95 


2.6 


2.5 


0.8 


Beamsize (natural weight) 


3.6" X 2.5" 


8.2" X 6.6" 


4.8" X 2.7" 


12.3" X 7.2" 


11.9" X 11.7" 


5.1" X 2.5" 


31.4" X 6.7" 


16.0" X 10.0" 


Beam position angle 


-6° 


1° 


7° 


8° 


-11° 


11° 


-3° 


18° 


Phase calibrator (B1950) 


0402-362 


0402-362 


1127-145 


1127-145 


0530+135 


2311-452 


2311-452 


0113-118 












& 2227-394 






Bandpass calibrator (B1950) 


0402-362 


0402-362 


1127-145 


1127-145 


0530+135 
& 3C 48 


2255-282 
& 3C 48 


2255-282 


0113-118 


Line and Continuum Properties 


















Peak continuum surface brightess (mjybeam" ) 


19.6 


45.8 


5.08 


9.13 


13.4 


13.6'' 


40.2 


27.3 


Total continuum flux density (mjy) 


110 


124 


123 


128 


20 


55 


55 


54 


Noise in continuum image (/ijybeam^^) 


35 


33 


25 


22 


20 


33 


45 


23 


Noise per image channel (/^Jybeam^^) 


120 


140 


66 


68 


56 


80 


135 


56 


Peak line surface brightness (mjybeam"^) 


< 0.36 


0.47 ±0.08 


< 0.20 


< 0.17 


< 0.24 


< 0.17 


Integrated line flux or 3(T to 6tT upper limit (Wm~ ) 


< 4 X 10-^^ 


(4.5 it 1) X 10-^^ 


< 0.47 X 10-^^ 


< 0.6 X 10-^^ 


< 0.76 X 10-^^ 


< 0.57 X 10-^^ 


Measured line width (FWHM; kms""') 




339 ± 75 














Deconvolved line width (FWHM; kms~^) 




215 ± 50 














Peak 8.3 GHz continuum (mjybeam"^) 


7.0 

















After offline Hanning smoothing. 
^Thc peak continuum strength for sources A & B is 11.5 mjybeam^^ and for C & D is 13.6 mjybeam""'^ in the CnB image. 
'^Assuming a FWHM equal to the velocity resolution of the data. 
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3. The Circinus Galaxy 

This remarkable low-latitude spiral is at a distance of 
about 4Mpc (thus 1" is 19 pc) and is the closest Seyfert 
galaxy (type 1) with many signs of a star burst in the nu- 
cleus. 

The H I distribution has been studied with the ATCA 

(Jones ct al. 1999) who find a total size of about 80' or 
90kpc with a total HI mass of 7 x IO^Mq. Detailed imag- 
ing of the circumstellar ring with a radius of about 220 pc 
in Ha (Elmouttic ct al. 1998a) shows a rotational speed of 
about 350kms~^ in rough agreement with an HI feature 
(Jones et al. 1999). This feature may lie close to the inner 
Lindblad resonance. 

A number of studies of the radio continuum have been 
carried out (Elmouttie et al. 1998b, Harnett et al. 1990, 
Elmouttic ct al. 1995). The peculiar radio lobes orthogo- 
nal to the disk of the galaxy are polarized and the nuclear 
source (about 0.1 Jy at 1.4 GHz) has a flat spectral index 
of q; = —0.06 {S oc i^") with a size limit of ^ 20 pc. Higher 
resolution radio images by Davies et al. (1998) at 8.4 GHz 
and 5 GHz with a resolution of 1" to 2" show a compact 
source with a size of about 0.6" or 12 pc. There is also a 
reported radio core of 19mJy at 2.3 GHz observed with a 
beam of 0.1" with the Parkes-Tidbinbilla Interferometer, 
which can be associated with the AGN. Davies ct al. sug- 
gest that the flat spectral index of the nucleus could arise 
from free-free absorption. 

There have been several studies of the molecular con- 
tent of the Circinus galaxy by Curran et al. (2001) and 
Elmouttie et al. (1998c). The former authors have carried 
out a systematic study of numerous molecules (eg four 
isotopes of CO in several transitions, CS, H2CO, HCN, 
HCO+ and other molecules). They propose that the H2 
density near the nucleus is in the range 2 x lO'^cm"^ to 
10^ cm~^. The molecular emission is distributed in a disk 
of radius 300 pc which is the likely source of the gas for the 
star formation. H2O masers have been extensively stud- 
ied by Greenhill et al. (2003a, 2003b). Doppler shifts up 
to ~ 460kms~^ are observed. 

In the optical and near IR, high resolution investi- 
gations of the Circinus galaxy have been carried out by 
Wilson et al. (2000) who have used HST imaging in [O HI], 
Ha, H2 and continuum bands in the optical and the near 
IR. A one-sided ionization cone is observed in Ha and 
[O III] and the nuclear starburst activity is shown in Ha 
images at distances of tens of parsccs from the centre. 
The circumstellar star- forming ring of radius ~ 200 pc is 
prominent. 

3.1. Results 

The ATCA continuum and line images, integrated spec- 
trum, and position-velocity diagram are shown in Figs. 1 
to 3. The measured line and continuum parameters are 
given in Table 1. 

The continuum image shows a single extended con- 
tinuum component. The continuum emission is predomi- 



nantly non-thermal, with a spectral index of -0.65 between 
1.4 GHz and 8.3 GHz (Elmouttie et al. 1998b). 

Line emission was detected in the nuclear region, co- 
incident with the peak continuum emission with a tail 
of emission extending 25" (500 pc) towards the north- 
east. The total area of line emission is 1.9 beam areas, or 
(250pc)2. 

The H91a -|- H92a spectrum integrated over the line- 
emitting region shows a clear line detection with 26 a sig- 
nificance with centroid at 480kms~^, compared to the HI 
systemic velocity of 449kms~^ (Juraszek et al. 2000). The 
line FWHM is 260kms~^ after deconvolving the instru- 
mental velocity resolution of 42kms^"'^. 

The position- velocity diagram (Fig 3) shows a gradient 
in the east-west direction of 130kms~^ over the central 
25" (475 pc) consistent with a disk rotating about a north- 
south axis with the east side receding. The position angle 
for the position- velocity slice was chosen by inspecting the 
first-moment image, which showed the velocity gradient to 
be east-west. Interestingly, the gradient direction is oppo- 
site that of the Br7 or H2 (Miiller Sanchez et al. 2006), the 
H2O masers (Greenhill et al. 2003) and the HI disk (Jones 
et al. 1999), all of which show the south-west side reced- 
ing. The gradient direction of the RRL emission is the 
same as that of the molecular outflow seen by Elmouttie 
et al. (1997) and studied by Curran et al. (1999). Thus, 
the H92a emitting gas is most likely associated with the 
outflow. 

The velocity span of the signal in the position- velocity 
diagram (130kms~^) is half as wide as the fine integrated 
over a 38" x 38" region in Fig 2 (260 km s"!). This indi- 
cates that there is extended low-level RRL emission out- 
side the region in the position-velocity diagram, i.e. at 
higher or lower declination. The line width integrated over 
the large area (260kms~^) is comparable to that of the 
molecular outflow (±190kms~^) and of the ionized out- 
flow (180kms-i; Veilleux & Bland-Hawthorn 1997) and 
is much larger than the width of the Br7 and H2 emission 
(30kms~^) seen on 8pc scale (Miiller Sanchez et al. 2006), 
further supporting the association of the RRL-emitting 
gas with the outflow. 

3.2. Modelling the Ionized Gas 

Two types of models have been considered for the RRL 
emission from the nuclei of external galaxies: one based 
on a uniform slab of ionized gas and the other based on 
a collection of compact H II regions. Such models have 
been discussed by Anantharamaiah et al. (2000) and ref- 
erences therein, and are documented in detail by Mohan 
(2002). These models take as constraints the integrated 
RRL strength at one or more frequencies, the radio con- 
tinuum spectrum and the geometry of the line emitting 
region. 

For modelling the RRL emission from the Circinus 
galaxy, we used the observed line strength (3.2 mJy), 
line width (260kms~^), size of the line-emitting region 



6 



A. L. Roy et al.: Recombination Line from Circinus Galaxy and NGC 1808 



line emission / mJy/beam * Icm/s 
20 40 60 




141316 



04 



12 10 08 06 
RIGHT ASCENSION (J2000) 

Peak continuum brightness = 98.5 mJy/beam 
Contours at -4, -2, -1,1,2,4, 8, ... 64 mJy/beam 

Fig. 1. ATCA 8.3 GHz + 8.6 GHz continuum image of tlie 
Circinus galaxy (contours) superimposed on the grey scale 
zeroth-moment image showing H91a + H92a line emis- 
sion. Peak continuum brightness is 98.5mJybeam~^ and 
contours are at -4,-2,-1,1,2,4,8,16,32, and 64mJybeam~^. 
Grey-scale peak = 65mJybeam~^ kms"-^. Beamsize 
is 11.0" X 9.4" at a P. A. of -76°, rms noise is 
0.12 mJy beam"-*^ channel"^. 



(250 pc), continuum emission (225 mJy), and spectral in- 
dex (-0.65) to constrain conditions in the ionized gas. 

We used the collection of H H regions model from 
Anantharamaiah et al. (1993). We considered a grid of 
models with electron temperature, Te, in the range 1000 K 
to 12 500 K, electron density, rie, in the range 10cm~^ to 
lO^cm"^, and effective diameter of the line-emitting gas 
(~volume"^/^), I, in the range 0.01 pc to 100 pc. 

We found that 3 x 10^ Mq to 1 x lO'^ Mq of ionized 
gas with Te ~ 5000 K, rie ~ 500 cm'^ to 5 x 10^ cm'^, 
and with a total effective diameter of the line-emitting 
gas of 3pc to 50 pc produced good matches to the line 
and continuum emission. We found that almost all of the 
RRL emission is due to stimulated emission amplifying 
the non-thermal continuum. The fraction of stimulated 
emission was around 90% in all models. Parameters de- 
rived for typical allowed models are given in Table 3. 
The allowed range of values is large since we are us- 
ing a single line strength measurement to constrain mul- 
tiple parameters. Much tighter constraints would come 
from a multi-transition study, as was done for Arp220 
by Anantharamaiah et al. (2000). 

The inferred mass of ionized gas and the Lyman con- 
tinuum flux required to maintain the ionization are sum- 
marized in Table 3. The flux is equivalent to the Lyman 



-3 
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in 
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■a 




-200 200 400 600 800 1000 1200 
Heliocentric velocity / km/s 

Fig. 2. ATCA H91a -I- H92a hnc profile integrated over 
the line-emitting region in the Circinus galaxy. Region of 
integration is a box of size 38" x 38" centred on RA 14 13 
10.01 dec -65 20 19.1. 
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Fig. 3. ATCA position- velocity diagram of H91q! + H92q: 
emission in the Circinus galaxy along a slice at con- 
stant declination (-65° 20' 19.1") through the maximum 
of the RRL emission, showing possible rotation. Beamsize 
is 11.0" X 9.4" at a P.A. of -76°, giving 10.6" resolution in 
RA, 36kms~^ in velocity (no Banning smoothing), and 
averaging over 9.8" in declination. 



continuum output of 300 to 9000 stars of type 05, which 
infers a star-formation rate of 0.2 yr~^ to 6Mq yr~^ 
when averaged over the lifetime of OB stars. 

This can be compared to star formation rates de- 
rived from other indicators following Hopkins et al. (2003). 
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Taking the peak 1.4 GHz continuum surface brightness of 
445 mJy beam" 1 in the 20" x 19" beam of the ATCA at 
the nucleus of Circinus (Elmouttic ct al. 1998b) yields 
a 1.4 GHz luminosity in the central 400 pc diameter of 
9.4 X 10^*^ Wm~^ and a corresponding star formation rate 
of 0.86 M© yr-i. The IRAS 60 fiia and 100 /im flux dcn- 
sitites yield a far-infrared (FIR) star formation rate of 
O.17M0 yr~^. The Ha-based star formation rate could not 
be estimated due to a lack of published Ha spectropho- 
tometry. The U-band magnitude of 12.87 (de Vaucouleurs 
et al. 1991) along with the Balmer decrement from Oliva et 
al. (1994) yields a U-band star formation rate of 1.63 M© 
yr~^. These estimates all agree well with the star forma- 
tion rate of 0.2 Mq yr~^ to 6 Mq yr~^ estimated from the 
RRL emission. 

4. NGC 1808 

NGC 1808, at a distance of 11 Mpc, is a spiral galaxy 
(Sbc pec) undergoing a star burst, with a FIR luminos- 
ity of 2 X 10^'^ Lq. This galaxy exhibits a bright nu- 
cleus, a bar, several optical hotspots (Sersic & Pastoriza 
1965), and dusty radial filaments over the central few kpc 
(Veron-Cetty & Veron 1985). The peculiar morphology 
and the starburst activity are attributed to tidal interac- 
tions with the companion NGC 1792 (Dahlem et al. 1990). 
HI observations show a strong outflow (Koribalski ct al. 
1993a), co-spatial with the dusty filaments. There is evi- 
dence both for and against Seyfert activity at the nucleus 
(Forbes, Boisson & Ward 1992; Kotilainen et al. 1996). CO 
and HI imaging (Dahlem et al. 1990; Koribalski, Dickey 
& Mebold 1993b; Koribalski et al. 1996) shows a cen- 
tral concentration of gas surrounded by a ^ 25" rotat- 
ing molecular ring, believed to be coincident with the 
inner Lindblad resonance of the kpc-scale bar. There is 
also evidence for non-circular gas motion (Saikia ct al. 
1990). High resolution radio continuum images (Saikia et 
al. 1990, Collison et al. 1994) show that the central kpc 
contains a population of compact sources. NIR contimuim 
and line imaging at 1.8" (Krabbc, Sternberg & Gcnzcl 
1994) and 1.0" (Kotilainen et al. 1996) shows a number of 
Br7, [Fell], and H2 knots. The Br7 and radio continuum 
knots are spatially correlated at a resolution of 1.8" but 
arc not exactly coincident at a higher resolution of 1.0". 
There is no correlation between the positions of the op- 
tical hotspots and the radio/Br7 knots, hence the former 
seem to trace regions of low dust extinction rather than 
regions of star formation. Most of the radio knots have a 
steep spectral index and are probably SNRs with a ther- 
mal fraction not exceeding ~ 30% at cm wavelengths. 
Using spectrophotometric models, Krabbe, Sternberg & 
Genzel (1994) and Kotilainen et al. (1996) derive an age 
of greater than 40 Myr to 50 Myr for the nuclear starbiu'st 
and much younger ages of less than 10 Myr to 15 Myr for 
the circumnuclear hotspots. Tacconi-Garman, Sternberg 
& Eckart (1996) imaged the galaxy in the NIR K band 
at 0.6" resolution and detected a number of unresolved 
sources, believed to be young massive super star clus- 



ters. Assuming a decaying star-formation rate, they de- 
rive ages similar to the earlier work for the nucleus and 
the hotspots. However, they also derive a much larger age 
of ~ 200 Myr for all of the sources assuming a history of 
continuous star formation. 



4.1. Results 

The brightest continuum source seen in Fig 4 is the nu- 
cleus, and the second brightest, labelled 'hotspot' was seen 
previously in the 1.4 GHz image of Saikia et al. (1990). 
Based on the 1.4 GHz and 8.3 GHz continuum flux densi- 
tites measured with a 3.1" beam, the spectral indices of 
the nucleus and the hotspot are calculated to be —0.57 
and —0.65 respectively. Higher resolution (sub-arcsecond) 
images show that the hotspot consists of two compact 
distinct sources. These arc referred to as AlO and All 
by Colhson et al. (1994). Though diffuse Br7 emission is 
seen in a region encompassing both sources, only a sin- 
gle Br7 knot is seen in this region and it lies very close 
to AlO (Krabbe, Sternberg & Genzel 1994, Kotilainen et 
al. 1996). Hence there seems to be more thermal gas to- 
wards AlO than All, even allowing for dust extinction. 
The radio continuum emission from AlO is weaker than 
from All in the high resolution images. However, images 
at a lower resolution of 3.1" at 8.3 GHz (this work) and 
1.4 GHz (Saikia et al. 1990) show that the smooth emis- 
sion peaks much nearer the position of AlO, indicating 
the presence of substantial diffuse emission around this 
source. 
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Fig. 4. 8.3 GHz continuum images of NGC 1808 with the 
VLA in CnB configuration. The synthesized beam is 3.6" x 
2.5" at a P. A. of —6°. The rms noise is 35 /iJybeam"-'^ and 
the contour levels plotted are (10, 20, 40, 60, 100, 120, 140, 
160, 180, 220, 260, 300, 340, 380, 420, 460, 500, 540)xrms. 
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H92a line emission integrated over the central 18" x 
20" was detected in the DnC configuration data (Fig 5), 
extending between 750 kms"'^ and 1100 kms"-'^. This ve- 
locity range is consistent with the HI velocity range of 
(800 to 1200)kms-i (Koribalski et al. 1993a), consider- 
ing the low signal-to-noisc ratio of the RRL detection. The 
observed FWHM of the line is (339 ± 75) kms'^ and the 
deconvolved line width, after correcting for ofHine Hanning 
smoothing and spectral resolution, is about 210 km s"-'^. 
The peak line emission (470 /uJybeam"^) is only 3.9 a. 
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Fig. 5. The Hanning smoothed H92a spectrum in NGC 

1808 for the DnC configuration data, integrated over a box 
of size of 18" X 20" centred on RA 05 05 58.58 dec -37 34 
36.5. The synthesized beam is 8.2" x 6.6" and the Icr noise 
in the channel image is 140 /iJy beam" ^. The peak of the 
line integrated over the region is 0.47 mJy and the contin- 
uum emission integrated over the same region is 101 mJy 
giving a line-to-continuum ratio of 0.005. 



4.2. Modelling the ionized gas 

We modelled the thermal gas in NGC 1808 using the mea- 
sured continuum and line strengths at 8.3 GHz and a 
spectral index of -0.7 at 8.3 GHz. We modelled the Hne- 
emitting gas as a collection of H II regions. For each combi- 
nation of input parameters in a given model and assuming 
conservatively that the thermal gas was spread homoge- 
neously over the entire region, we corrected the intensity 
of the non-thermal radiation for absorption by the ther- 
mal gas and calculated the expected free-free continuum 
emission from the model H II regions. The resulting free- 
free emission was constrained to be less than the observed 
thermal emission of 30 mJy, which was estimated from the 



total observed continuum emission of 101 mJy assuming a 
thermal fraction of 30 %. 

The allowed range of electron density is 100 cm""^ to 
1000 cm~'^ and the effective diameter of the line-emitting 
gas (~volume^/^), I, is in the range 7pc to 300 pc. The 
ionization rate is 3x10^^ s~^ to lxlO'''^s^^. The free-free 
flux density at 8.3 GHz is < 1.0 mJy to 20mJy which is 
< 1 % to 20 % of the total continuum. 

Model solutions were calculated for Te of 5000 K, 7500 
K, and 10 000 K. The derived parameters varied by much 
less than the range in the allowed parameter values for 
any given Tc value. Thus, we present parameters for Tc = 
5000 K, but the values for = 7500 K or 10 000 K would 
be almost the same. The fraction of stimulated emission, 
as for the Circinus galaxy, was around 90 % in all models. 

Using the derived value of iVLyc= 3 x 10^^ s~^ 
to IxlO^^s^^, the Br7 flux is expected to be 
3 x 10~^^ Wm"2 to 9 X 10~^^Wm~^. For comparison, 
Krabbe, Sternberg & Genzel (1994) measured the Br7 flux 
from the hotspot region to be 2.6xl0~^^ Wm~^ and the 
data of Kotilainen et al. (1996) yield 1.6xlO-i'^Wm-2. 
Both measurements are within the broad range expected 
given the (extinction-free) RRL strength. 

To estimate the effect of small changes in the 6„ and 
(3 coefficients on the derived physical parameters of the 
ionized gas, which is a concern since the line emission 
is strongly affected by stimulated emission, we note that 
when we ran models for adjacent lines having roughly 
the same observed parameters, we obtained essentially the 
same values of derived parameters and conclude that the 
effects of changes are small. 

Parameters derived for typical allowed models are 
given in Table 3. 

5. The Non-Detections 

The H92q! line was not detected in the other galaxies ob- 
served with the ATCA (NGC 6221, IC5063, IRAS 18325- 
5926, and NGC 7552) or with the VLA (NGC 4038/9, 
IIZw40, NGC 7552, and VV114). Upper limits were 
0.8 mJy beam" 1 to 1.4 mJy beam" ^ (3cr) for the ATCA 
and 0.17mJy beam"^ to 0.24mJy beam"^ (3ct) for the 
VLA observations (Tables 1 and 2). The continuum im- 
ages are shown in Fig 6. 

IC 5063 is classed as a NLRG, having 5 x 10^3 W Hz^i 
at 1.4 GHz. The nuclear component was resolved into a 4"- 
long linear triple structure at 8.6 GHz by Morganti et al. 
(1998). Our 8.4 GHz image (Fig 6a) with the larger beam 
of 6.4" X 5.8" did not resolve the jet, as expected. 

We used the upper limits on the H91a -I- H92a 
line emission to constrain the properties of the nuclear 
emission-line regions. We assumed that the line FWHM 
was the same as the velocity resolution of the data (from 
Tables 1 & 2). We modelled each source as a collection 
of uniform-density spherical HII regions. The model so- 
lutions were constrained to produce an 8.3 GHz free-free 
emission equal to 30 % of the continuum flux density inte- 
grated over the central region and an H92q! line emission 
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equal to the observed upper limit. The derived value for 
the gas density varied between 100 cm"^ and 5 x 10^ cm~^ 
and the total effective diameter of the H II regions varied 
between 5pc to 200 pc with the effective size decreasing 
with increasing density. The total flux of ionizing photons 
required to maintain such nebulae on the limit of our de- 
tectability ranged between 3 x 10^^ s~^ and 2 x 10^'*s~^ 
depending on density and size of the individual HII re- 
gions. This provides upper limits to the number of 05 stars 
powering the ionization in the range 100 to 70 000, with 
the higher values corresponding to those models whose 
combination of parameters yield high free-free continuum 
emission from the ionized gas. 

Parameters derived for typical allowed models are 
given in Table 3. 

NGC 6221 contains two emission-line regions in the 
nucleus - a low-excitation H II region probably photoion- 
ized by hot. young stars and a faint high-excitation re- 
gion with relatively broad emission lines characteristic of 
Seyfert 2 nuclei (Veron et al. 1981). Our continuum image 
at 8.3 GHz (Fig 6b) shows emission from the nucleus, bar 
and spiral arms, consistent with that seen at 1.4 GHz by 
Koribalski & Dickey (2004). 

We used the upper limits on the H91a -|- H92a 

line emission to constrain the properties of the nuclear 
emission-line regions following the approach described for 
IG 5063. Parameters derived for typical allowed models are 
given in Table 3. 

IRAS 18325-5926 is an ultraluminous infrared 
galaxy with a compact central radio component (size 
< 0.1") which accounts for one fifth of the total emis- 
sion at 2.3 GHz (Roy et al. 1994). Our 8.3 GHz image in 
Fig 6c is consistent with that of Bransford et al. (1998), 
showing an unresolved component with no significant ex- 
tended structure. The flux density in our observation, 
30 mJy ± 4 mJy, is lower than the 36.8 mJy ±0.3 mJy mea- 
sured by Bransford et al. (1998) (after interpolating be- 
tween their 4.8 GHz and 8.6 GHz measurements). The 
source might therefore have varied between 1993 Jul and 
1995 Nov. 

We used the upper limits on the H91q; -|- H92a 
line emission to constrain the properties of the nuclear 
emission-line regions following the approach described for 
IC 5063. Parameters derived for typical allowed models are 
given in Table 3. 

NGC 4038/9 is a pair of interacting galaxies with ex- 
tended tidal tails and numerous star clusters, including su- 
per star clusters (Whitmore & Schweizer 1995, Whitmore 
et al. 1999). Neff & Ulvestad (2000) imaged the multi- 
frequency radio continuum emission in this galaxy at high 
resolution (1.0" to 2.6") and identified a host of compact 
sources, some of which exhibit a spectral index flatter than 
—0.4, indicating a high thermal gas fraction. Our contin- 
uum image (Fig 6d) is consistent with theirs, showing the 
same collection of compact sources immersed in diffuse 
emission. Our image is at lower resolution and has a fac- 
tor two higher noise. 



We used the upper limit on the RRL emission to 
constrain the properties of the strongest thermal source 
54.96- 06.1 (sec Fig 6d). Neff & Ulvestad (2000) deter- 
mined the diameter of this source to be 70 pc (0.7") and 
derived an rms density of ~ 400 cm~^. 

We modelled this source as a collection of uniform- 
density spherical H II regions, constrained by the H92a 
upper limit, the 8.3 GHz flux density of Neff & Ulvestad 
(2000) as an upper limit, the observed 4.8 GHz to 8.3 GHz 
spectral index, and the infra-red excess (IRE), LfirZ-^Lyc 
(derived from Lfir=1.3x1036 W; Neff & Ulvestad (2000) 
and the A^Lyc of the models) to exceed unity. No solutions 
were obtained for the thermal fraction of the 8.3 GHz con- 
tinuum /th > 0.6 if the RRL constraint was included. The 
derived values for the gas density and the size of the H II 
region varied between 5xl0^cm~'^ to lO'^cm"^ and 40 pc 
to 5pc, and the IRE was 6 to 10. The H92a line strength 
is '--^ 80 % of the 3cr upper limit if the thermal fraction is 
not less than 30% and hence might be detectable with 
sensitive observations. However, the observed flat spec- 
tral index of —0.26 ± 0.13 cannot be reconciled with the 
low derived thermal fraction unless there is additional low 
density thermal gas which does not emit 8.3 GHz RRL. 
It is surprising that we do not detect RRL emission from 
other regions in the galaxy. Since it is unlikely that the 
starburst has declined rapidly over the last million years, 
a search for RRLs at much lower resolution, possibly using 
the GBT, might prove more fruitful. 

IIZw40 was imaged by Beck et al. (2002) at 15 GHz 
using the VLA and discovered four unresolved objects. 
They identified the continuum emission as free-free and 
suggested that these are similar to the young 'supernebu- 
lae' discovered in NGC 5253 (Turner, Beck & Ho 2000) and 
He 2-10 (Kobulnicky & Johnson 1999). They show that 
these supernebulae account for almost all of the 12 /zm 
emission from the inner core. We show our continuum im- 
age in Fig 6e, which, as expected due to the lower reso- 
lution, shows the multiple sources as a single unresolved 
component. 

We used the upper limits on the H92q: + H93a line 
emission with a model consisting of an individual spherical 
H II region of uniform density, constrained to produce a 
15 GHz continuum free-free emission of 0.6 mJy. Solutions 
were obtained for rie > 5000 cm~^. The expected 8.3 GHz 
H92a line flux for the various model solutions for all four 
objects is 60 % to 70 % of the observed 3a upper limit 
for densities between 5000 cm""^ and lOOOOcm^'', is about 
40 % for rie ~ 50 000 cm~'^ and decreases further for higher 
densities. 

NGC 7552 contains three star forming rings of radii 
1.0 kpc, 1.9 kpc, and 3.4 kpc at the co-rotation and 
Lindblad resonance radii of the bar. The inner-most ring is 
seen in radio continuum (Forbes, Kotilainen & Moorwood 
1994). Our continuum image (Fig 6f) shows the ring struc- 
ture seen previously by Forbes et al. A number of Br7 
knots lie along the inner ring which show weak correla- 
tion with radio knots, but these are not spatially resolved 
in our observations. 
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We used the upper limits on the H91a + H92a 
line emission to constrain the properties of the nuclear 
emission-line regions. Parameters derived for typical al- 
lowed models are given in Table 3. 

VV 114 is undergoing a merger and the two nuclei 
VV114E and VV114W arc separated by 15" (5.6 kpc). 
The E and W components are barely resolved in the con- 
tinuum image shown in Fig 6g. 

The model for the thermal component of VV 114 E was 
constrained by the upper limit to the H92a line emission 
and a free-free emission of fth x Sc where Sc is 27.2 mJy 
at 8.3 GHz. The IRE was constrained to be greater than 
unity by assuming that the FIR luminosity of VV 114 E is 
1x10^^ Lq. Consistent solutions were obtained only for a 
thermal fraction < 20%, i.e., Sc < 5.4 mJy. If this thermal 
gas is optically thin at 8.3 GHz, then the ionization rate 
is < 3x10^^8-^ 



6. Conclusions 

We have detected HQla or H92a lines in emission in the 
Circinus galaxy and NGC 1808 with integrated flux densi- 
ties of 3.2 mJy and 0.47 mJy respectively using the ATCA 
and the VLA. We established upper limits on the RRL 
strength in NGC 4038/9, IIZw40, NGC 6221, NGC 7552, 
IRAS 18325-5926, IC 5063, and VV 114. The detected line 
strengths infer ionized gas masses of 3000 Mq to 10^ Mq 
in the Circinus galaxy and 10^ in NGC 1808, and cor- 
responding star formation rates of 0.2 Mq yr~^ to 6Mq 
yT~^ and 0.3 yr"-'^ to 0.6 Mq yr"-'^, depending on the 
model conditions. The star formation rate estimated from 
the RRL detection in Circinus agrees well with rates es- 
timated from radio, FIR, and U-band luminosities using 
previously-calibrated relations. 

Since the detectability of RRL emission depends so 
sensitively on the gas density, future prospects for using 
RRL emission as diagnostics of star formation activity 
are good if one observes multiple transitions and includes 
higher frequencies. Multi-transition studies are important 
to provide sensitivity to a wide range of gas densities and 
the details of which particular transitions are detected 
provide the density distribution of the ionized gas. The 
line strengths have been predicted and found to increase 
with frequency, making studies at 22 GHz, 43 GHz, and 
86 GHz attractive due to their larger line-to-continuum ra- 
tios. Further, the interpretation of line strengths at higher 
frequency is simplified since stimulated emission is much 
less important towards high frequency, giving a direct re- 
lationship between RRL emission strength and the ionized 
gas mass, and hence to the star formation rate. 

Acknowledgements. The National Radio Astronomy 
Observatory is a facility of the National Science Foundation 
operated under cooperative agreement by Associated 
Universities, Inc. The Australia Telescope Compact Array 
is part of the Australia Telescope, which is funded by the 
Commonwealth of Australia for operation as a National 
Facility managed by CSIRO. 



1^ 



H :2 > H H H H 

P r ° » ^ o 
g "< S S o n 

So" - a 5 ^ 



N 2. &- S 

<t S 2 B 



A A A2>S g 
8§x S 3 X 



^ 
o X „ 

Ol |_L Q 







A ^ A Ol Ol ~a 
P P g ^ 



o c c 



o 3 « 



JO SJ o c 



a ft > 



2 a 
9 " 

o o 



a a a cfi 



B B B B 
S' S' 5' S' 



lA A A S E S E 



3 3 3 B 



3 3 3 3 
« Cfi « K 



H 5? > H H H H 
3- o p o »; ^ ^ 

3 o" - a ft s 

§ 0- <■ § g 

<—' ui p m 

5' E5. g- 

o a 0) » 



CO I-' W W Oi CJ1 

O X X X o o 





CI 

c 


1 


1,0 


pc 


to 


° 




Ol 

o 


X ^ 




o 
o 
o 






n 


^ =^ 

























A -vj CO rf:^ I 

-2 X ^ 

- 5 ° ^ ° 



or Ol 
O O 
. O O 



sag 



A A A 5 g 



O 

o 
o 



A. L. Roy et al.: Recombination Line from Circinus Galaxy and NGC 1808 



11 



References 

Anantharamaiah, K. R., & Goss, W. M. 1997, RMxAC, 6, 58 
Anantharamaiah, K. R., Viallefond, F., Mohan, N. R., Goss, 

W. M., & Zhao, J. H. 2000, ApJ, 537, 613 
Anantharamaiah, K. R., Zhao, J. H., Goss, W. M., & 

ViaUefond, F. 1993, ApJ, 419, 585 
Baars J. W. M., Genzel, R., PauHny-Toth, I. I. K., Witzel, A. 

1977, A&A, 61, 99 
Beck, S. C., Turner, J. L., Langland-Shula, L. E., Meier, D. S., 

Crosthwaite, L. P., Gorjian, V. 2002, AJ, 124, 2516 
Bell, M. B., & Seaquist, E. R. 1978, ApJ, 223, 378 
Bell, M. B., Seaquist, E. R., Mebold, U., Reif, K., & Shaver, 

P. A. 1984, A&A, 130, 1 
Bransford, M. A., Appleton, P. N., Heisler, C. A., Norris, R. 

P. & Marston, A. P. 1998, ApJ, 497, 133 
Churchwell, E., & Shaver, P. A. 1979, A&A, 77, 316 
CoUison, P. M., Saikia, D. J., Pedlar, A., Axon, D. J., Unger, 

S. W. 1994, MNRAS, 268, 203 
Cornwell, T. J., Uson, J. M. & Haddad, N. 1992, A&A, 258, 

583 

Curran, S. J., Rydbeck, G., Johansson, L. E. B., & Booth, R. 

S. 1999, A&A, 344, 767 
Curran, S. J., Johansson, L. E. B., Bergman, P., Heikkila, A., 

& Alto, S. 2001, A&A, 367, 457 
Dahlem, M., Aaho, S., Kleiu. U.. Booth, R., Mebold, U., 

Wielebinski, R. & Lesch, H. 1990, A&A, 240, 237 
Davies, R. E., Forbes, D. A., Ryder, S., Ashley, M. C. B., 

Burton, M., Storey, J. W. V., Allen, L. E., Ward, M. J., & 

Norris, R. P. 1998, MNRAS, 293, 189 
de Vaucouleurs et al. 1991, Third Reference Catalogue of 

Bright Galaxies, version 3.9 
Elmouttie, M., Haynes, R. F., Jones, K. L., Ehle, M., Beck, R., 

& Wielebinski, R. 1995, MNRAS, 275, L53 
Elmouttie, M., Haynes, R. F., & Jones, K. L., 1997, PASA, 14, 

140 

Elmouttie, M., Koribalski, B., Gordon, S., Taylor, K., 

Houghton, S., Lavezzi, T., Hayues, R. F., & Jones, K. L. 

1998a, MNRAS, 297, 49 
Elmouttie, M., Haynes, R. F., Jones, K. L., Sadler, E. M., Ehle, 

M. 1998b, MNRAS, 297, 1202 
Elmouttie, M., Krause, M., Haynes, R. F., & Jones, K. L. 

1998c, MNRAS, 300, 1119 
Forbes, D. A., Boisson, C, & Ward, M. J. 1992, MNRAS, 259, 

293 

Forbes, D. A., Kotilainen, J. K., & Moorwood, A. F. M. 1994a, 
ApJ, 433, L13 

Forbes, D. A., Norris, R. P., Williger, G. M., & Smith, R. C. 

1994b, AJ, 107, 984 
Frater, R. H., Brooks, J. W. & Whiteoak, J. B. 1992, JEEEA, 

12, 103 

Greenhill, L. J., Kondratko, P. T., Lovell, J. E. J., Kuiper, T. 
B. H., Moran, J. M., Jauncey, D. L., & Baines, G. P. 2003a, 
ApJ, 582, Lll 

Greenhill, L. J., Booth, R. S., EUingsen, S. P., Herrnstein, J. 
R., Jauncey, D. L., McCulloch, P. M., Moran, J. M., Norris, 
R. P., Reynolds, J. E., & Tzioumis, A. K. 2003b, ApJ, 590, 
162 

Harnett, J. I., Whiteoak, J. B., Gardner, F. F., Reynolds, J. 

E., & Tziousmis, A. 1990, MNRAS, 244, 130 
Hibbard, J. E., van der Hulst, J. M., Barnes, J. E., & Rich, R. 

M. 2001, AJ. 122. 2969 



Hopkins, A. M., Miller, C. J., Nichol, R. C, Connolly, A. J., 
Bernardi, M., Gomez, P. L., Goto, T., Tremonti, C. A., 
Brinkmann, J., Ivezic, Z., & Lamb, D. Q. 2003, ApJ, 599, 
971 

Joues, K. L., Koribalski, B. S., Elmouttie, M., & Haynes, R. 

F. 1999, MNRAS, 302, 649 
Juraszek, S. J., Staveley-Smith, L., Kraan-Korteweg, R. C, 

Green, A. J., Ekers, R. D., Haynes, R. F., Henning, P. A., 

Kesteven, M. J., Koribalski, B., Price, R. M., Sadler, E. 

M., & Schroder, A. 2000, AJ, 119, 1627 
Koribalski, B., Dahlem, M., Mebold, U., & Brinks, E. 1993a, 

A&A, 268, 14 

Koribalski, B., Dickey, J. M. & Mebold, U. 1993, ApJ, 402, 
L41 

Koribalski, B., Dettmar, R. -J., Mebold, U., & Wielebinski, R. 

1996, A&A, 315, 71 
Koribalski, B. & Dickey, J. M. 2004, MNRAS, 348, 1255 
Kobulnicky, H. A., & Johnson , K. E. 1999, ApJ, 527, 154 

(KJ99) 

Kotilainen, J. K., Forbes, D. A., Moorwood, A. F. M., van der 

Werf, P. P., & Ward, M. J. 1996, A&A, 313, 771 
Krabbe, A., Sternberg, A., & Genzel, R. 1994, ApJ, 425, 72 
Mohan, Niruj R., Anantharamaiah, K. R. & Goss, W. M. 2001, 
ApJ, 557, 659 

Mohan, Niruj R. 2002, PhD Thesis, IISc, Bangalore, India 
Morganti, R., Oosterloo, T., Tsvetanov, Z. 1998, AJ, 115, 915 
Miiller Sanchez, F. Davies, R. I., Eisenhauer, F., Tacconi, L. 

J., Genzel, R., & Sternberg, A. 2006, A&A, 454, 481 
Neir, S. G. & Ulvestad, J. S. 2000, AJ, 120, 670 
Oliva, E., Salvati, M., Moorwood, A. F. M., & Marconi, A. 

1994, A&A, 288, 457 
Phookun, B., Anantharamaiah, K. R., & Goss, W. M. 1998, 

MNRAS, 295, 156 
Rodn'guez-Rico, C. A., Goss, W. M., Zhao, J.-H., Gomez, Y. 

& Anantharamaiah, K. R. 2006, ApJ, 644, 914 
Roy, A. L., Norris, R. P., Kesteven, M. J., Troup, E. R. & 

Reynolds, J. E. 1994, ApJ, 432, 496 
Roy, A. L. 1995, PhD Thesis, University of Sydney 
Roy, A. L., Goss, W. M., Mohan, Niruj R., Anantharamaiah, 

K. R. 2005a, A&A, 435, 831 
Roy, A. L., Goss, W. M., Mohan, Nhuj R., Oosterloo, T. 

& Anantharamaiah, K.R. 2005b, in The Evolution of 

Starbursts. 331st Wilhelm and Else Heraeus Seminar. 

16 - 20 August, 2004, Bad Honnef, Germany, Susanne 

Huettemeister, Eva Manthey, Dominik Bomans, Kerstin 

Weis (eds) AIP Conf Proc 783 (AIP, New York), p303 
Saikia, D. J., Unger, S. W., Pedlar, A., Yates, G. J., Axon, 

D. J., Wolstencroft, R. D., Taylor, K. and Gyldenkerne, K. 

1990, MNRAS, 245, 397 
Seaquist, E. R., & Bell, M. B. 1977, A&A, 60, LI 
Sersic, J. L., & Pastoriza, M. 1965, PASP, 77, 287 
Shaver, P. A. 1978, A&A, 68, 97 

Shaver, P. A., Churchwell, E., & Rots, A. H. 1977, A&A, 55, 

435 

Tacconi-Garman, L. E., Sternberg, A., & Eckart, A. 1996, AJ, 
112, 918 

Turner, J. L., Beck, S. C, & Ho, P. T. P. 2000, ApJ, 532, L109 
(TBHOO) 

Vcron-Cetty, M. P., & Vcron, P. 1985, A&A, 145, 425 
Veron, M. P., Veron, P., & Zuiderwijk, E. J. 1981, A&A, 98, 
34 

Veilleux, S., & Bland-Hawthorn, J. 1997, ApJ, 479, L105 
Whiteoak, J. B. 1970, ApL, 5, 29 



12 



A. L. Roy et al.: Recombination Line from Circinus Galaxy and NGC 1808 



Whitmorc, B. C, & Schwcizcr, F., 1995, AJ, 109, 960 
Whitmore, B. C, Zhang, Q., Leitherer, C, Fall, S. M., 

Schweizer, F., & Miller, B. W. 1999, AJ, 118, 1551 
Wilson, A. S., Shopbell, P. L., Simpson, C, Storchi-Bergmann, 

T., Barbosa, F. K. B., & Ward, M. J. 2000, AJ, 120, 1325 
Wright, A. E. 1974, MNRAS, 167, 273 
Wright, A. E. & Otrupcek R. 1990, Parkes Catalog, ATNF 
Zhao, J. -H., Anantharamaiah, K. R., Goss, W. M. & 

Viallefond, F. 1996, ApJ, 472, 54 



A. L. Roy et al.: Recombination Line from Circinus Galaxy and NGC 1808 



13 




Fig. 6. (a) 8.4 GHz continuum image of IC 5063 using the ATCA. The synthesized beam is 6.4" x 5.8" at a P.A. of 24°. 
The peak brightness is ISOmJybeam"^ and the contour levels are at -1.8, 1.8, 3.6, 7.2, 14.4, 28.8 and 57.5mJybeam~^. 
(b) 8.4 GHz continuum image of NGC 6221 using the ATCA. The synthesized beam is 11.8" x 7.9" at a P.A. of 32°. 
The peak brightness is 20.2mJy beam^i and the contour levels are at -0.4, 0.4, 0.8, 1.2, 1.6, 2.0, 2.4, 2.8, 4.8, 6.9, 8.9, 
10.9, 14.9, 14.9, 17.0 and 19.0mJybeam-^ (c) 8.3 GHz continuum image of IRAS 18325-5926 using the ATCA. The 
synthesized beam is 6.3" x 5.6" at a P.A. of 50°. The peak brightness is 30.9m.Jy bcam"-'^ and the contour levels are 
at -0.3, 0.3, 0.6, 1.2, 2.5, 4.9, 9.9 and 19.8mJy beam"!, (d) 8.3 GHz VLA continuum image of NGC 4038/9 using the 
C configuration. The rms noise in the image is 25 /uJybeam"^ and the contour levels are (5, 10, 15, 20, 25, 30, 35, 40, 
50, 60, 70, 80, 100, 120, 140, 160, 180, 200)xrms in image. The synthesized beam is 4.8" x 2.7" at a P.A. of 7°. (e) 
8.3 GHz continuum image of 11 Zw 40 using the VLA in the D configuration. The synthesized beam is 11.9" x 11.7" 
at a P.A. of —11° and the rms in the image is 20;uJybeam~^. The contour levels are (5, 7, 10, 15, 20, 40, 60, 80, 100, 
150, 200, 250, 300, 350, 400, 450, 500, 550, 600, 650) times the rms. (f ) 8.3 GHz contimmm image of NGC 7552 using 
the VLA in the CnB configuration with uniform weighting. The rms in the image is 58/xJybeam~^ and the contour 
levels are (10, 20, 30, 40, 50, 60, 70, 80, 90, 100, 110, 120, 130)xrms. The synthesized beam is 3.1" x 1.7" at a P.A. 
of 20°. (g) 8.3 GHz continuum image of VV 114 in the D configuration. The contour levels are (10, 20, 40, 80, 160, 
320, 500, 700, 900, 1100, 1170)x rms in the image, which is 23AtJy beam~i. The synthesized beam is 16.0" x 10.0" at 
a P.A. of 18°. 



